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The finite-size spectrum for the multi-channel Kondo model is derived analytically from the 
exact solution, by mapping the nontrivial Z„ part of the Kondo scattering into that for the RSOS 
model coupled with the impurity. The analysis is performed for the case of n — 2S = 1, where n is 
the number of channel and S is the impurity spin. The result obtained is in accordance with the 
Kac-Moody fusion hypothesis proposed by Affleck and Ludwig. 
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The multi-channel Kondo problem is one of the most 
interesting issues in strongly correlated electron systems 
It possesses the non- trivial fixed point character- 
ized by a non-Fermi liquid behavior for the overscreening 
case. A possibility of its realization in some Uranium 
compounds has been discussed extensively ||||- Re- 
cently, Affleck and Ludwig developed a boundary confor- 
mal field theory (CFT) approach to this problem, and de- 
rived various non-trivial results for low-energy properties 
j7j. Their key proposal is the so-called fusion hypothesis; 
i.e. when absorbing the impurity spin into conduction 
electrons, the selection rule for spin quantum numbers 
is changed following the Kac-Moody fusion rules. Using 
this hypothesis they predicted the finite-size spectrum, 
which plausibly yields the ir/2 phase shift for n < 25, as 
well as the non-Fermi liquid spectrum for n > 25, where 
n is the number of channel (flavor) , and 5 is the impurity 
spin. 

On the other hand, the Bethe-ansatz solutions to var- 
ious Kondo models have been obtained [§J|, which in 
principle can be compared with the above prediction. 
This was performed for the single channel Kondo model, 
confirming that the Bethe-ansatz indeed gives the spec- 
trum in accordance with the fusion hypothesis |]l0f . For 
the overscreening model however, the correct finite- 
size spectrum has not been obtained analytically from the 
Bethe-ansatz solution. The difficulty comes from that 
the multi-channel model relies on the string hypothesis 
even for the ground state, which is valid only in the ther- 
modynamic limit |^|,^| . Hence the ordinary Bethe-ansatz 
approach suffers from this pathology, leaving the analytic 
computation of the finite-size spectrum still open. 

In this paper, we propose a way to deduce the finite- 
size spectrum of the overscreening Kondo model analyt- 
ically from the exact solution. A key idea is to con- 
vert the nontrivial Z„ part of the original Kondo scat- 
tering into that for the critical restricted solid-on-solid 
(RSOS) model coupled with the impurity. The finite- 
size spectrum is then calculated by the functional equa- 
tion method based upon fusion hierarchy plldT^]. The 



obtained spectrum plausibly reproduces the fusion hy- 
pothesis of Affleck and Ludwig. 

We start by recapitulating the difficulty in the ordi- 
nary Bethe-ansatz approach. The multi-channel Kondo 
model is given by 



H = H e kc{ la C kla + J 



k,k' ,l 7 a,<j' 



J kla 



(CW • S)c k 'la' (1) 



with the antifcrromagnetic coupling J > 0, where I is 
the flavor index (I — 1, 2, • • • , n), a =f, j is the spin in- 
dex, and other notations are standard. This model was 
solved by the coordinate Bethe-ansatz method, which has 
described the thermodynamics exactly However, a 

naive application of finite-size techniques ]l4| , |l5| based on 
the string hypothesis turns out to fail for the overscreen- 
ing case, and gives only the gaussian part of the spectrum 
for the spin sector. Namely, the total spectrum leads to 
a pathological form, 
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for the finite system with linear size N, where v is the 
Fermi velocity, and uq, n s , and nf are non- negative 
integers. Here the first term expresses the charge ex- 
citation specified by the quantum number Q, and the 
second is the spin excitation labeled by the magneti- 
zation S z . Since the spin sector of the multi-channel 
Kondo model may be correctly described by the level- 
n SU(2) Wess-Zumino-Witten model with the central 
charge cy/zw — 3rt/(2 + n) Q, it is seen from (^) 
that the Z n parafermion sector with central charge c = 
cwzw — 1 = 2(n — l)/(n + 2) is lacking. This type of dif- 
ficulty is common to many integrable models with affine 
symmetry Eg ]. 

We thus resort to alternative techniques other than the 
coordinate Bethe ansatz to recover the correct spectrum 
for the lacking Z„ parafermions. Though it seems not 
easy to resolve this problem completely, we propose an 
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analytic way to investigate the nontrivial Z n parafermion 
part. For this purpose, we first exploit the following prop- 
erties of the S-matrix for "physical particles" [l7|-[Tg|] in 
the overscreening Kondo model. The bulk S-matrix in 
the spin sector is decomposed into two parts 0,^9) ; the 
S-matrix for the 6-vertex model which is equivalent to 
the gaussian theory, and that for the Z n model. Since the 
Z„ model can be described in terms of the RSOS model 
the corresponding S-matrix is given by the face 
weight in the RSOS model |l^,^2|. A remarkable point 
for the overscreening model is that the interaction be- 
tween "physical particles" and the impurity is described 
by the S-matrix for multi-kinks, which is given by the 
fusion of the face weights of the RSOS model with the 
fusion level p — n — 2S. This has been previously noticed 
by Fendley and Martins [p3[ . Therefore, it is seen 
that nontrivial properties in the overscreening model is 
essentially determined by the RSOS model coupled with 
the impurity. Based on these observations, we now con- 
sider the transfer matrix T(u) for the RSOS model with 
the impurity, 

T(u) = W(u\{a j+1 })...W(u\{a N })W p (u\{a N+1 }) 

xW(u\{ai})...W(u\{*j-i}), (3) 

which could reproduce the nontrivial properties of the 
multichannel Kondo model. Here V^'s are the face 
weights of the RSOS model with a spectral parameter 
u, and {a} is the set of spins around the face. The im- 
purity counterpart of the face weight W p is given by the 

fusion of w ^MMMM, 



W p (n\{a}) = J2U W 

{<j'}k=l 



p+l 



k A 



W , (4) 



where A = n/(n + 2). Note that from this equation 
one can indeed reproduce the Z n sector of the thermody- 
namic Bethe-ansatz equation obtained for the overscreen- 
ing case |J|. 

We wish to study the above RSOS model with the im- 
purity to deduce the correct spectrum of the overscreen- 
ing model, and show how the fusion hypothesis emerges 
in this framework. With a view to avoiding a pathology 
in the string hypothesis, we deal with (||) following the 
functional equation method pj--|l3|] . The essence of this 
method is that one can seek for the solution to the eigen- 
value problem by observing analyticity properties of the 
transfer matrix. We shall slightly generalize it to take 
into account the effect of the impurity. Let us first con- 
struct the fusion hierarchies of this model following the 
standard fusion procedure [^6|,^5). Introduce first the 
fused face weight for host electrons |^5],[l2 24 , 

W^(u\{a}) = [] s^(u) J2f[w(u+(k- l)A|{er'}), 



k=0 



fc=l 



(5) 



MI ^"sinA ' an d the fused face weight for 



where Sk(u) 

the impurity, which is given by eq. (|5|) with Sk (tt) replaced 

by 8 k (u) = ITj=o sin ( u +( k + E 2 1 - j) X )/ sin A - We then 
obtain the transfer matrix for the fusion hierarchies, 



T q (u) = W lq ...W 1 ' q W p ' q W lq ...W 1 ' q . 



(6) 



From the generalized Yang-Baxter equations satisfied by 
W 1,q and W p ' q , it is easily seen that the commutativ- 
ity of the transfer matrix holds, [T q (u),T q («)] = 0, 
which guarantees the construction of the fusion hierar- 
chies. The fusion hierarchies satisfy the following func- 
tional equations p4 12|, 



TqTI = /_!/! + / T 2 
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where T q = T(u + gA), and 



u = 



siiu // + </A i ^ N j-j sin(w + (q 

3=0 



sin A 



2 



sin A 



(7) 

(8) 



(9) 



Following Kliimper and Pearce |L2j, we recast the above 
functional equations into the following form, 



tg*? = (l + t? _ )(l + *o )» l<?<n-l (10) 

with t q = T q ^T q+1 1 f-if q , and t% = = 0. If we omit 
the face weight for the impurity, W p , the analyticity strip 
of t q {u) is given by @ 
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2(n + 2) 



2(n + 2) 



(11) 



Note that t q (u) has poles in this strip. For a large 
system with linear size N, t q (u) can be described by 
t q {u) = (z q (u)) y q (u) , where (z q (u)) N expresses the 
host electron part and y q (u) the impurity part. From 
eq.(|l(]), these functions respectively satisfy 



z q (u)z q (u + A) 
z q - 1 {u + \)z<i+ l {u) 



(12) 



for 1 < q < n — 1, with the condition z°(u) = z n (u) = 1, 
and 



V ° V \i = 1, Kg<n-1, 



(13) 



with y°=yft = l. 

Now our task is to find the solution to z q (u) and as well 
as to y q (u), which should respect the correct pole struc- 
ture of t q (u) in the analyticity strip (pi]). The solution 
to eq.(12) was already found [Ejfl, 
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z q {u) 



sin[s±%)] 

a[ n±2 {u + qX) y 



(14) 



which indeed possesses the pole structure correctly. This 
immediately yields the well-known finite-size spectrum of 
Z n parafermion theory Q , 
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where m = — 2j, — 2j + 2 , . . . , 2 j — 2,2j, and j = 
0,s,l, Note that only the holomorphic piece of 

conformal weights appears since the Kondo problem is 
effectively described by a chiral system . 

In order to obtain the solution to the impurity part, 
y q (u), we have to know how the presence of the impurity 
changes analyticity properties of t q . Although this prob- 
lem may not be easy in general, it is found that for the 
special case of p = 1 , the analyticity strip ( pi] ) is slightly 
modified due to the impurity. Hence the unique solution 
to y q can be easily obtained, which reproduces the pole 
structure of t q under the condition j/q = yfi = 1. We find 
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In the following we will restrict our discussions to the 
case p = 1, owing to the above technical reason. 

Let us now turn to the calculation of the finite-size 
spectrum. To this end, it is convenient to introduce the 
new functions, a q (x) = [fi^ix- (q + f )J0] _1 - Apply- 
ing the method of Kliimper and Pearce p2], we obtain 
the non- linear integral equation for a q , 



In a 9 = lne q + J2 K Q.r*H 1 



D q , 



(17) 



where the asterisk stands for the convolution. Here D q 
is the vector independent of x (see eq.(p6|)) and Fourier 
components of the matrix K are given by the solution of 
the equation, 

(Ko(z)5 m ,l — <5m,i-l — Sm,l+l)Kl, n = —5m,n+l — ^m,n-\i 

(18) 

with Kq(z) = — 2 cosh(wz/n). lne 9 is determined by 
t q = (z q (u)) y q (u) , of which the bulk part is given in 
ref. [Lq ] . The impurity part is given by 

from which one can easily obtain the scaling limit of In e 9 , 

lim lne q {±(x + In TV)) = (bulk part) =p (20) 
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We are now ready to compute the finite-size spectrum 
from the finite-size corrections to the eigenvalue of the 
transfer matrix, 
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where T2 .. is the finite-size corrections to the transfer 
matrix. The integral equation for b q (x) is easily found 

IB3, 



m& 9 



^ir 9 , r *ln(l + a r ) + C 9 , 



(22) 



where the kernel if is given by 



sin ^ sin — 



1*1 



(23) 



and C 9 is the vector consisting of constant elements 
which do not contain 1/N corrections. From eqs.(|l7|) 
and fl22), one can obtain the finite-size spectrum 12 1, 
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q=l ' q=l 

where L(x) is the Rogers dilogarithm. Here aj^. 
lim^^oo a q (x), and = lim^^l + a q (x)), which 
evaluated from the asymptotics of eq. ( ]l0| ) as , 

sin 2 (9 .„ sin 2 (g + l)6» 



sing^ sin(<7 + 2)9' 



A q 



svaqO sin(<7 + 2)9 



(25) 



where 9 is a multiple of n/(n + 2). D q is determined from 
the asymptotic behavior of cq. (|l7|) , 



D q = (bulk part) - 



2pqni 



(26) 



The second term comes from the impurity contribution, 
and is responsible for determining the selection rule for 
quantum numbers, as we will see shortly. 

We thus find from eqs.(|24]) and ( |2^ ) the correction of 
the finite-size spectrum due to the impurity, 
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Here s — 1,2, ...,n + 1, and m = s — 1 (mod 2), 
\m\ < s — 1 = 2j. The second term in the bracket can be 
evaluated from the braid limit p2] , p7[ , 
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sin A 



sin(u — A/2) 
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T(u) = 2 cos ■ 



(28) 



Since the left-hand side of eq.(p8|) is evaluated as inm/n, 
the second term in the bracket of eq.(p7f) is equal to iirm. 
Consequently we end up with the simple formula, 



2tt 

^E R sos = ir T m P- 



(29) 



This completes the calculation of the finite-size correc- 
tions due to the impurity. 

Now combining the ordinary finite-size spectrum for 
Z n parafermions eq.(|l^) and the impurity contribution 
eq.(B9), we obtain the desired result, 



NErsos + (m+p) 2 
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n + 2 



in 



const, 



(30) 



where m = 2j (mod2), and j = 0, 1/2, 1, n/2. One 
can see that the spectrum fits in with Z n parafermion 
theory, and only the selection rule for quantum numbers 
is changed via m — > (m + p), after incorporating the im- 
purity effect. We think that this may be a microscopic 
description of the fusion hypothesis. 

To confront our results with those of Affleck and Lud- 
wig 0| , let us look at the total finite-size spectrum, which 
is given by the sum of cqs.(||) and ( ^o|). One can easily 
see that the term — (m + p) 2 /in in (|30|) should be can- 
celed by S 2 /n in (|J) by suitably choosing the quantum 
number for S z . This in turn modifies the quantum num- 
bers in SU(2) n Kac-Moody algebra. As a result, the total 
spectrum is reduced to 



NE = (Q - n) 2 | j(j + 1) 
2nv An n + 2 

+n Q + n s + rif, 



(flavor part) 



(31) 



with newly introduced quantum numbers j = \j —p/2\ 
where j = 0, 1/2, 1, n/2. Note that Q and j are the 
charge and spin quantum numbers for free electrons with- 
out Kondo impurities. Consequently we can say that the 
sole effect of the Kondo impurity is to modify the selec- 
tion rule for quantum numbers in the spin spectrum, via 
j — ► j, which indeed leads to non- Fermi liquid properties 
for the overscreening Kondo effect. This is the essence 
of the Kac-Moody fusion hypothesis proposed by Affleck 
and Ludwig 0. 

In this paper we have been concerned with the spe- 
cial case p = (n — 25) = 1 in the overscreening model. 
As mentioned above, we have encountered some technical 
problems for treating the case with general p, which have 
not been resolved yet. Nevertheless we think that the so- 
lution to y q for general p is given by eq. (|l6|) , and the 
corresponding finite-size spectrum should take the form 
of eq. (|3l|) . We hope to clarify these points in the future 



work to establish the present conclusion for the general 
case. 
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